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Tau neutrino and antineutrino interactions with nucleons in large underground or under ice detec- 
tors will be important signals of astrophysical and atmospheric sources of neutrinos. We present here 
a theoretical update of the deep inelastic scattering contribution to the tau neutrino and antineu- 
trino charged current cross sections with isoscalar nucleon targets and proton targets for incident 
neutrinos and antineutrinos in the energy range from 10 GeV to 10 TeV. Next-to-leading order 
quantum chromodynamic corrections, target mass corrections and heavy quark effects are included. 
Uncertainties in the cross section associated with the structure functions a low momentum transfers, 
the input parton distribution functions, scale dependence and flavor number scheme are discussed. 
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I. INTRODUCTION 



Neutrinos are signals of hadronic interactions that pro- 
duce mesons, followed by meson decay. There are many 
models which predict neutrino production in astrophysi- 
cal sources [l[ as well as in the atmosphere Q. IceCube 
[3| and the detector element DeepCore |4( and large un- 
derwater detectors [5] aim to measure neutrino signals 
from atmospheric production, from individual sources 
and the isotropic neutrino flux summed over all the 
sources. 

The phenomenon of neutrino oscillations is an impor- 
tant ingredient in evaluating the potential signals in de- 
tectors. From combined solar and atmospheric neutrino 
measurements, we know that the three flavors oscillate 
into each other 6]. Direct production of tau neutri- 
nos in the atmosphere is small [7|, but oscillations of 
atmospheric muon neutrinos to tau neutrinos can result 
in large upward v T + v T fluxes at certain energies ja, |9j. 
Astrophysical distances are much larger than terrestrial 
distances, so even if the dominant production in astro- 
physical sources is v^ + v^, a substantial fraction of the 
flux will oscillate to v T + v r . With the low atmospheric 
backgrounds at high energies, tau neutrino signals from 
convential and exotic [lfj astrophysical sources may be 
important. The DeepCore detector may be capable of 
distinguishing tau neutrino events |3]. 

An essential input to theoretical predictions and exper- 
imental analyses of tau neutrino induced events is the tau 
neutrino charged current cross section [llj . We provide 
here an update of earlier work [12L |l3[ with a summary 
of tau neutrino and antineutrino deep inelastic scatter- 
ing (DIS) charged-current cross sections in the energy 
range of 10 GeV - 10 TeV. Tau mass corrections, rel- 
ative to the muon neutrino charged current cross sec- 
tion, are important for some of this energy range. The 
low Q 2 extrapolations of the nucleon structure functions, 
not included in Refs. [12i LL3j, are also relevant to the 
low end of the energy range. We evaluate the next-to- 
leading order quantum chromodynamic corrections [14j . 
the effect of target mass corrections [15I - KL7J and the range 



of predictions from different ways of incorporating the 
charm quark contribution |18l420l |. A new estimate of 
the theoretical errors for the DIS cross sections for tau 
neutrino and antineutrino charged-current interactions is 
presented. 



II. CROSS SECTION FOR TAU NEUTRINOS 



In the lower end of the energy r ang e discussed here, 
targ et mass [15|-[l7J], tau mass fill [l2| and quark mass 
|l4| corrections are important |13| |. The low momentum 
transfer behavior of the structure functions is also a con- 
sideration. To discuss these corrections, we define the 
integration variables x = Q 2 /2p • q and y = p ■ q/p ■ k, 
where the momentum assignments are 



v/v T {k) 



N(p) -> r/f(k') 
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X 
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The target nucleon mass M comes in primarily through 
the replacement of the parton light cone momentum 
fraction x by the Nachtman variable 77 = x(2/(l + 
v /l + 4M 2 a; 2 /Q 2 )). Target mass corrections (TMC) also 
modify the structure functions themselves, where the un- 
corrected functions are denoted by i^ and the target mass 
corrected structure functions are F™ c . The difference 
between Fi and i^. TMC is largest at large x. For a sum- 
mary of target mass corrections to the structure func- 
tions, see, for example, Refs. |15l - fl7l ]. 

The tau mass m T is kept explicitly in the differen- 
tial cross section. By keeping terms proportional to 
m 2 /ME^, five of the six target mass corrected struc- 
ture functions allowed by the general Lorentz structure 
of the hadronic current appear in the differential cross 
section [Tljj . The expression for the tau neutrino charged 
current differential cross section, for deep inelastic scat- 



tering (DIS), is 



conveniently evaluated using 
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where neutrino (antincutrino) scattering requires 

+F TMC („ F TMC) incq ^ 

For muon neutrino charged current interactions, F^- MC 
and F^ MC are generally omitted since they are sup- 
pressed by a factor of at least m 2 /2ME v relative to the 
contributions of F\, F2 and F3. Albright and Jarlskog 
noted that at leading order, in the massless quark and 
massless target limit, F4 = and F5 = F2/2X [ll(. These 
are called the Albright- Jarlskog (AJ) relations. Correc- 
tions to these relations at next-to- leading order (NLO) 
in quantum chromodynamics including charm mass cor- 
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rections appear, for example, in Ref. 

A lepton mass correction to the charged current cross 
section comes from the kinematic limits on x and y [111 . 
Il2j |. For completeness, we reproduce them here: 



2M{E V ~ m T ) 
a — b < y < a 



< x < 1 
6, 
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where the quantities a and b are 
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A final kinematic issue relates to the DIS nature of the 
scattering. For quasielastic scattering, e.g., v T n — > rp, 
the structure functions are proportional to the delta func- 
tion S(W 2 — M 2 ) where W 2 is the invariant mass of the 
hadronic final state. These multiply the nucleon form fac- 
tors [2l|, |22j ■ Nuclear physics models are used to evaluate 
few pion production [2J, |24j . To avoid double counting 
these exclusive contributions in our DIS evaluation, we 
require that the invariant mass of the hadronic final state 
W be larger than a minimum value. Our standard choice 
is W m in = 1.4 GcV. In terms of x and Q 2 , this means 
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The minimum W 2 as well as the tau mass come into 
play in the differential distribution in y. This is most 
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where the limits of integration are 
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Q m . n = 2E 2 (l-e)(l-y)-m 2 T , 

= 2MEy + M 2 - W^ in = 2MEy - A 2 , 
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where e = \ 1 — m 2 / 1 (1 — y)E) 
Now y ranges between 



&ndA 2 = W 2 :-M 2 . 
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Even with W 2 > W^ in , Q 2 can be below the nom- 
inal minimum Q 2 for the parton distribution functions 
(PDFs) used to evaluate F 4 TMC . We use a phenomcno- 
logical parameterization of the structure functions first 
discussed by Capella, Kaidalov, Merino and Tran Than 
(CKMT) [5| and later used in Ref. 26] for the low Q 2 
extrapolations needed in v^N scattering. The CKMT 
parameterization of the structure functions at low Q 2 ac- 
counts for the nonperturbative behavior of i^™ . The 
CKMT parameterization is matched at Q 2 — Q 2 to the 
perturbative evaluations of i^ TMC with PDFs. Below this 
cutoff momentum transfer Q c , we use the A J relations 
for F^- MC and F^ MG , regardless of the incident neutrino 
energy. In our results below, we exhibit the sensitivity 
of the cross section to different choices of Q c . These 
results are quite similar to those using the Bodek- Yang- 
Park prescription [27j_for the low Q 2 extrapolation of the 
structure functions 26]. 

While we use the A J relations for low Q 2 , we can test 
the A J relations for higher Q 2 > Q 2 values. Setting 
Fj MC = for all Q 2 results in a neutrino cross section 
only 1% higher than for F 4 TMC ^ at 10 GeV, while 
for antineutrinos, the enhancement is 4%. At 100 GeV, 
the Fj MC contribution is negligible. By substituting 
f tmc = p™ G /2x at 100 GeV, the cross sections change 
by less than 1%. They differ by 2% for v r N and by 9% 
for v T N at 10 GeV as compared to the full calculation of 

pTMC 
^5 

TMC of Ref. 



in Ref. 16 1 . We use the full NLO calculation with 
llH for E v < 100 GeV and Q 2 > Q 2 C . 



The AJ relations for E v > 100 GeV reproduce the full 
calculations of the cross sections. 



III. RESULTS 

For the cross sections evaluated here we have used 
two different sets of parton distribution functions. The 



Energy [GeV] 


a vN [1(T 38 cm 2 ] 


gun [10~ 38 cm 2 ] 


10 


0.916 (1.26, 0.690) 


0.291 (0.574, 0.160) 


10 1.25 


3.77 (4.22, 3.44) 


1.48 (1.90, 1.21) 


10 1.5 


10.4 (10.9, 9.97) 


4.43 (4.93, 4.05) 


10 1.75 


23.7 (24.3, 23.2) 


10.6 (11.2, 10.2) 


10 2 


48.9 (49.5, 48.4) 


22.8 (23.4, 22.3) 


10 3 


5.69 x 10 2 


3.02 x 10 2 


10 4 


4.30 x 10 3 


2.76 x 10 3 



TABLE I: The i/ r -isoscalar nucleon and i^-isoscalar nucleon 
charged current cross section using the CKMT parameteriza- 
tion for Q 2 < 2 GeV 2 in the structure functions, matched to 
the perturbative structure functions with GJRF PDFs [2^ . 
using the charm mass corrections from Ref. [Ig] and target 
mass corrections from Ref. [la. Hq). Here, the factorization 
scale is set to /i — Q, and the hadronic invariant mass is 
limited to W > 1.4 GeV (M, 1.7 GeV) in the cross section. 



Energy [GeV] 


(Ju V [10~ 38 cm 2 ] 


ap p [10" 38 cm 2 ] 


10 


0.472 (0.646, 0.356) 


0.428 (0.840, 0.230) 


10 1.25 


2.12 (2.36, 1.94) 


2.08 (2.69, 1.69) 


10 1.5 


6.21 (6.49, 5.98) 


6.04 (6.78, 5.51) 


10 1 ' 75 


14.8 (15.1, 14.6) 


14.2 (15.0, 13.5) 


10 2 


31.6 (32.0, 31.4) 


29.8 (30.6, 29.1) 


10 3 


3.99 x 10 2 


3.72 x 10 2 


10 4 


3.28 x 10 3 


3.26 x 10 3 



TABLE II: The «/ T -proton and z? T -proton, as in Table I, for 
the DIS charged current cross section. The hadronic invariant 
mass is limited to W > 1.4 GeV (M, 1.7 GeV) in the cross 
section. 



Gluck, Jimenez-Delgado and Reya (GJR) PDFs 28] for 
a fixed flavor number scheme with three light flavors is 
used. One feature of the GJR PDFs is that the fits are 
done for Q 2 > 0.5 GeV 2 . This offers the option for a 
range of cut-off scales Q 2 , below which the phenomeno- 
logical CKMT parameterizations are used [25|, [26| . Our 
standard choice is Q 2 = 2 GeV 2 . We also use the 
CTEQ6.6M PDFs ;29]~. This set includes 44 variations 
of the best fit parameter values to help characterize the 
error associated with the PDF fit. 



The charged current cross sections divided by inci- 
dent tau neutrino or antineutrino energy with Q 2 = 2 
GeV 2 and the three flavor GJR PDFs with the factor- 
ization scale in the structure functions set to /j, 2 = Q 2 , 
are plotted in Fig. 1 for iso-scalar nucleon targets. The 
solid lines show <jcc{ v tN)/E and the dashed lines show 
&cc(v T N)/E. The cross sections for a few energies 
are also listed in Table 1. The upper curves show the 
cross sections when W m i n = M, while the lower curves 
have W m in = 1.4 GeV. Unless otherwise noted, our de- 
fault evaluation is with W m [ n = 1.4 GeV. Table 1 also 
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FIG. 1: The tau neutrino- (solid) and antineutrino- (dashed) 
isoscalar nucleon charged current cross section, for W > 1.4 
GeV (lower curves) and W > M (upper curves). Target mass 
corrections, NLO QCD and low Q extrapolations of the struc- 
ture functions below Q 2 = 2 GeV 2 are included in the evalu- 
ation, as in Table 1. 



lists results when W m in = 1.7 GeV. For completeness, 
we include in Table 2 the tau neutrino-proton and tau 
antineutrino-proton DIS charged-current cross sections. 
The rise in the cross section divided by incident neu- 
trino energy at low incident energy is largely from the 
kinematic threshold effect. The additional structure 
function F^f MC , which comes in with a minus sign, has 
a numerical effect even at 100 GeV incident energy. Ne- 
glecting F^f MC leads to an error of 13% for tau neutrinos 
and 29% for tau antineutrinos for E u — 100 GeV, with 
even larger errors at lower energies. At higher energy, W 
boson propagator already has an effect. For example, the 
neutrino cross section decreases by 6% for E = 10 3 GeV 
and by 30% at E = 10 4 GeV relative to the evaluation 
with a simple four-Fermi interaction. 



Overall, the kinematic suppression of tau production 
and the structure function term in the differential cross 
section lead to tau neutrino and antineutrino cross sec- 
tions lower than the muon neutrino cross section over a 
surprisingly wide range of energies. With the combined 
corrections from the structure functions and the kine- 
matic limits, the tau neutrino and antineutrino charged 
current cross sections are less than the corresponding 
muon neutrino cross sections by about 25% at 100 GeV 
incident neutrino energy. Even at 1 TeV, the tau neutrino 
(antineutrino) charged current cross sections are reduced 
by 5%(7%) compared to the muon neutrino and antineu- 
trino cross sections for an isoscalar target. We show in 
Fig. f5]the ratio of the charged current DIS cross sections 
for v T N and v T N scattering for lf m i n = 1.4 GeV. 

For scattering with a proton target and Wmin = 1.7 




0.8 



E [GeV] 

FIG. 2: The ratio of acc{v T N)/crcc{v^N) neutrino- (solid) 
and antineutrino- (dashed) DIS cross sections, for W > 1.4 
GeV. 



GeV, the neutrino ratio is between 4-7% lower than for 
what is shown in Fig. Hfor E=10-100 GeV. The antineu- 
trino ratio for a proton target with W m [ n = 1.7 GeV is 
slightly larger than the ratio in Fig. [2j At E=100 GeV, 
it is larger by about 2% and by ~ 1% at 10 3 GeV. 



A further demonstration of the kinematic suppression 
due to tau lepton production is seen in the ^-distribution. 
In Fig. [3l we show da/dy/E for tau neutrino and an- 
tineutrino scattering as well as for muon neutrino and 
antineutrino scattering for an isoscalar target. We show 
the differential distribution for two incident energies: for 
E = 10 GeV (upper figure) and E = 50 GeV (lower fig- 
ure). With muon neutrino scattering, one sees a nearly 
constant cross section as a function of y, because at 
these energies, the cross section is valence dominated, 
so da/dy ~ q(x,Q 2 ) for the valence parton distribu- 
tion function q(x,Q 2 ). At low y, the figure illustrates 
the consequence of the Wmin choice. For tau neutrino 
scattering, the limits in Q 2 and y cut off the high y 
distribution. For antineutrino scattering, valence dom- 
ination in the structure functions yields schematically 
da/dy ~ (1 — y) 2 q(x, Q 2 )- This means that the antineu- 
trino differential cross section is already falling with in- 
creasing y, so the kinematic effects of producing a tau 
lepton are less noticible. For both muon antineutrino 
and tau antineutrino scattering, the low-Q 2 range is more 
emphasized, since the average y value is lower than for 
neutrino scattering. 

The theoretical error on the DIS cross section has a 
number of components which we discuss in the remainder 
of this section. We begin with the low-Q 2 extrapolation 
using the CKMT structure functions. The GJR PDFs 
are fit to Qo = 0-5 GeV 2 . This allows us to compare the 
cross sections with choices for the transition point be- 
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(a) E= 10 GeV 




FIG. 3: The differential cross section da/dy/E for DIS 
charged current scattering with an isoscalar nucleon target, 
with W > 1.4 GeV for (a) 10 GeV and (b) 50 GeV incident 
energies. 



tween PDFs and the CKMT structure functions Q 2 . We 
find that for neutrinos, setting Q 2 — 0.5 GeV 2 results 
in <jcc{v t N) to less than 1% of the cross sections with 
Q c = 2 GeV 2 and W m [ n — 1.4 GeV. As noted, antineu- 
trino scattering occurs at lower values of Q 2 , so the value 
of Q 2 has somewhat more importance for antineutrinos. 
At E = 10 GeV, a(D T N) is about 3% lower for Q 2 C = 0.5 
GeV 2 than for Q 2 = 2 GeV 2 with W min = 1.4 GeV. By 
E = 100 GeV, the discrepancy is reduced to ~ 1%. 

The factorization scale dependence of the cross section 
introduces a larger theoretical error. For v T N scattering, 
setting the factorization scale to fi 2 = 0.5Q 2 enhances 
the cross section by close to 3%, while for fi 2 = 2Q 2 
decreases the cross section by approximately 4% as com- 
pared to /j 2 = Q 2 at E — 10 GeV. For D T N scattering, 
the range at 10 GeV is between ±5% The factorization 
scale dependence declines with energy. 



The choice of parton distribution function also has an 



IV. FINAL REMARKS 
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FIG. 4: An estimate of the theoretical errors for v T (triangles) 
and v T (squares) deep-inelastic charged current cross sections 
for selected incident (anti-)neutrino energies with W m i n — 1.4 
GeV. 



implication for the predicted DIS charged current cross 
section. Using the CTEQ6.6M PDFs in the evaluation 
of the cross section using Ref. [16] yields the v T N [p T N) 
cross section between 0.6-2.7% (0.9-3.0%) higher than 
using the three-flavor GJR PDFs for E = 10 - 10 4 GeV. 

The CTEQ6.6 PDFs are available with a set of 44 vari- 
ations of the best fit parameter values. Using these 44 
sets, we find that the resulting error in the cross section 
is largest for the lowest energy. For v T N scattering, the 
PDF uncertainty is 2.7% at 10 GeV, reducing to 1.4% at 
10 4 GeV. Slightly larger for v r N, the estimated CTEQ6.6 
PDF error is 3.6% at 10 GeV and 2% at 10 4 GeV. 

The final uncertainty in the DIS theoretical predic- 
tion considered here is due to the choice of calcula- 
tional scheme to include the charm quark. At E = 10 
GeV, charm contributions are negligible, but as the en- 
ergy increases, charm contributions become more impor- 
tant. Beginning at E — 100 GeV, we have evaluated 
ccci^rN) and acc{v T N) using the ACOT and the S- 
ACOT-x schemes [1{|, which we compare to the fixed 
flavor scheme result of Ref. [16| including massive charm 
quark corrections. 

Using the same CTEQ6.6M PDFs, we see a small range 
of theoretical predictions for the flavor schemes we used 
here. The S-ACOT- X result for v T N at E = 100 GeV is 
lower by close to 4% than the calculation using Ref. [l6| , 
and close to 6% lower for D T N. The discrepancy reduces 
to close to zero at E = 10 3 GeV, and it remains small at 
E = 10 4 GeV. 

Other schemes to incorporate charm are possible. A 
detailed discussion and comparison of these schemes ap- 
pears in Ref. 20]. The ultrahigh energy implications 
of the scheme dependence in neutrino-nucleon scattering 
are discussed in Ref. (30| . 



We have presented an evaluation of the tau neutrino 
and tau antineutrino cross sections for DIS charged- 
current scattering with isoscalar nucleon targets. The 
Albright- Jarlskog approximations for F^ and F$ are ex- 
cellent at high energies, however, at E = 10 GeV, for 
antineutrino scattering, the deviation from the full cal- 
culation of F$ MC is on the order of ~ 10% for tau an- 
tineutrino scattering. We have considered the numerical 
implications of a range of factorization scales, of the PDF 
choice and of the cutoff scale Q 2 C below which a low-Q ex- 
trapolation is used for the structure functions. We have 
also made an estimate of the uncertainty of the flavor 
scheme dependence of the inclusion of the charm quark 
contribution. 

A summary of our estimates of the theoretical uncer- 
tainties to the DIS cross section is shown in Fig. 21 where 
the various contributions have been added in quadrature. 
At 10 GeV, the largest error comes from the scale de- 
pendence of the structure functions, while at the highest 
energies, a more important element is the choice of PDF. 
At intermediate energies, the flavor scheme (S-ACOT-% 
versus the fixed flavor scheme) for the incorporation of 
charm gives the largest variation in the theoretical pre- 
diction. 

Our evaluation of the DIS cross section relies on a min- 
imum value of the final state hadronic invariant mass 
W . The inclusion of few pion and quasi-elastic neutrino 
scattering is required for a full evaluation of the neutrino- 
nucleon cross section. At 10 GeV, the quasi-elastic cross 
section is about 1/3 of the neutrino DIS cross section 
with W > 1.4 GeV, and comparable to the antineutrino 
cross section for the same W cutoff. The few pion cross 
sections are lower. Combining the DIS, exclusive and 
quasi-elastic contributions in v^ scattering appears in, 
e.g., Ref. [3l|. This procedure is applicable to tau neu- 
trino scattering as well. 

The tau neutrino and antineutrino cross sections have 
been measured by the DONuT Collaboration, with re- 
sults summarized in Ref. [321 ] . For an average of tau 
neutrino and antineutrino cross sections at an average 
energy of 115 GeV, 

a avg /E = 0.39 ± 0.13 ± 0.13 x 10~ 38 cm 2 GeV _1 



for the 9 v T + D T identified interactions. The OPERA 
experiment has likely seen its first tau neutrino induced 
event from the CERN neutrino beam [33|]. While the 
DONuT measurement is the first direct measurement of 
the tau neutrino cross section, one expects that tau neu- 
trino charged current interactions will play an important 
role in detailed measurements of neutrino oscillations in 
the Earth, for example, with the DeepCore detector in 
IceCube 0. 
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